Stability of adhesion clusters under constant force 



T. Erdmann 1 and U. S. Schwarz ,2 

1 Max- Planck- Institute of Colloids and Interfaces, Potsdam, Germany 

^Institute of Theoretical Physics, University of Leipzig, D-04103 Leipzig, Germany 

We solve the stochastic equations for a cluster of parallel bonds with shared constant loading, 
rebinding and the completely dissociated state as an absorbing boundary. In the small force regime, 
cluster lifetime grows only logarithmically with bond number for weak rebinding, but exponentially 
for strong rebinding. Therefore rebinding is essential to ensure physiological lifetimes. The number 
of bonds decays exponentially with time for most cases, but in the intermediate force regime, a 
small increase in loading can lead to much faster decay. This effect might be used by cell-matrix 
adhesions to induce signaling events through cytoskeletal loading. 



Biological systems have to be able to change quickly 
in response to external stimuli. This is one of the rea- 
sons why molecular bonds in biological system are based 
on non-covalent interactions and have short lifetimes of 
the order of seconds. In order to achieve long-lived as- 
semblies, cells in multicellular organisms adhere to the 
extracellular matrix and to each other through clusters 
of adhesion molecules. The number of receptors in adhe- 
sion contacts can range from just a few (e.g. for tethering 
of leukocytes to vessel walls or in the nascent contacts 
close to the leading edge of a locomoting cell) to ~ 10 5 
(e.g. in mature cell-matrix contacts). Most types of ad- 
hesion clusters are coupled to the cytoskeleton and have 
to function under mechanical load, which leads to expo- 
nentially increased dissociation rates Rebinding of 
broken bonds is often facilitated by the densely packed 
arrangement of molecular bonds in the cluster. How- 
ever, if the cluster operates under force, it is likely to 
be pulled away by some elastic relaxation process in the 
moment the last bond has been broken. Therefore the 
cluster usually cannot rebind from the completely disso- 
ciated state. Recently, the interest in the role of force at 
adhesion clusters has strongly increased also in the bi- 
ological community, since it has been shown that force 
at cell-matrix adhesions correlates with contact size and 
intracellular signaling 

Quantitative characterization of adhesion bonds has 
made tremendous progress during the last decade, mainly 
on the level of single molecules 0, El Due to the 
low binding energies of adhesion bonds, thermal activa- 
tion is important and theoretical models are required to 
interpret experimental data |(|. In order to make con- 
tact with situations of biological interest, the quantita- 
tive effort now has to be extended to clusters of adhesion 
bonds. Clusters also allow to study the effect of rebind- 
ing, which is difficult to address on the level of single 
molecules jjj. Recently, micropipette techniques have 
been used to study cluster dissociation under a linear 
ramp of force ||, in good agreement with a theoretical 
analysis by Seifert ||. However, physiological loading of 
adhesion clusters is usually more or less constant on the 
timescale of cluster lifetime. The stability of adhesion 
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FIG. 1: Schematic representation of an adhesion cluster un- 
der constant force: there are No = 5 receptor-ligand pairs, 
of which i = 3 are closed and equally share the constant 
dimensionless force /. Single closed bonds rupture with dis- 
sociation rate k — koe^ % and single open bonds rebind with 
force- independent association rate k on . Our model has three 
parameters: cluster size No, dimensionless rebinding rate 
7 = k on /ko and dimensionless force /. 

clusters under constant force has been first modeled by 
Bell 0, but his treatment was based on a simplifying 
deterministic equation for the mean number of bonds. 

In this Letter, we present exact and simulation results 
for a stochastic version of the Bell-model. In contrast 
to the deterministic model, the stochastic one allows to 
treat the completely dissociated state as an absorbing 
boundary. Moreover, it includes fluctuation and non- 
linear effects which are important for small adhesion clus- 
ters. The main objective of this work is to obtain analyt- 
ical results for generic features of adhesion clusters. We 
present several new formulae for cluster lifetime as a func- 
tion of cluster size, rebinding rate and force, which now 
can be used for quantitative analysis of adhesion experi- 
ments. Although we do not address the specific features 
of cell-matrix adhesions, our model suggests an appealing 
mechanism by which cells might induce signaling events 
through cytoskeletal loading. 

Following Bell JlJ , we consider a cluster with a constant 
number Nq of parallel bonds. At any given time, each of 
the different bonds can be either open or closed. The 
constant force F applied to the cluster is assumed to be 
shared equally between the i closed bonds (0 < i < N ). 
We assume that a single bond under force F ruptures 
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with the dissociation rate k = koe F ^' lFb introduced by 
Bell |l|, which can be rationalised by modelling bond 
rupture as thermally activated escape over a sharp tran- 
sition state barrier |6(. In this framework, the force scale 
Fb = kBT/xb is set by thermal energy fc^T and the dis- 
tance Xb between the potential minimum and the tran- 
sition state barrier along the reaction coordinate of rup- 
ture. For typical values xj, ~ 1 nm and T ~ 300 K, we 
find the typical force scale Fb ~ 4 pN. For the single 
bond association rate k oni we assume that it is indepen- 
dent of force. Physiological values for both ko and k on 
are expected to be in the 1/s-range. For the following, 
it is useful to introduce dimensionless time r = hot, di- 
mensionless rebinding rate 7 = k on /ko and dimensionless 
overall force / = F/Fb- Since bond rupture is a discrete 
process, the stochastic dynamics of the bond cluster can 
be described by a one-step Master equation ^(J 

dpi 

—± =r (i + l)p i+ i + g(i - l)p t -i - [r(i) + g(i)]pi (1) 

where Pi(r) is the probability that i bonds are closed 
at time r. The reverse and forward rates between the 
possible states i are 



r(i) = ie f/ \ g(i) = -f(N - i) 



(2) 



Depending on the experimental setup, rebinding from 
the completely dissociated state (i = 0) might be pos- 
sible (reflecting boundary) or not (absorbing boundary). 
For f — and a reflecting boundary at i = 0, we deal 
with natural boundaries, that is given reasonable initial 
conditions, no special equations are needed to treat the 
boundaries. For finite /, r(0) = is required to prevent 
i from becoming negative. An absorbing boundary at 
i = requires g(0) — 0. The mean number of closed 
bonds N as a function of time r is N = (i) = YlitoiPi- 
From the Master equation Eq. QJ, one can derive 



^> =-<r(i)> + te(i)> 



(3) 



If r and g are both linear functions in i, Eq. © becomes 
an ordinary differential equation for TV. This suggests to 
study the following deterministic equation 



^- = -NeV N + 7 (N -N) 
dr 



(4) 



as has been done by Bell 1] for constant loading and by 
Seifert Q for linear loading. However, for finite force / 
solution of Eq. (J2J does not give the correct result for the 
first moment, since then the rate r defined in Eq. |J2j is 
non-linear in i and the average in Eq. © cannot be taken. 
More importantly, a differential equation like Eq. (@J fol- 
lows from the stochastic equations only in the case of 
natural boundaries. In order to treat the biologically rel- 
evant case of an absorbing boundary, one therefore has 
to study the stochastic description Eq. QJ. 
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FIG. 2: Solution of the Master equation gives the state prob- 
abilities pi for i bonds being closed at time r (0 < i < No). 
Here they are plotted for iVo = 10. (a) 7 = and / = 0. 
(b) 7 = 1 and / = for a reflecting boundary at i = 0. (c) 
Same for absorbing boundary, (d) 7 = and / = 50. (a) 
and (b) follow from Eq. JjjJ, (c) is obtained from Monte Carlo 
simulations, and (d) follows from Eq. ©. 



A full stochastic solution amounts to finding the set 
of state probabilities pt (r) as a function of the three di- 
mensionless parameters No, 7 and /. For force / = 
and a reflecting boundary, the solution results from the 
generating function given by McQuarrie [ll| : 

p . (T) - (N \ h + e-^Y(l-e-^) N "- 1 (5) 

(l + 7 )iVo • W 

Here and in the following we use the initial condition 
iV(0) = Nq. If also rebinding 7 = 0, then we deal with 
the simple case of independently decaying bonds. Fig. [5^, 
shows that in this cluster with 10 bonds decays 

from % = 10 to by visiting each of the intermediate 
states to an appreciable degree. In order to stabilize the 
cluster, one has to introduce rebinding. Then there is 
fast relaxation to a stable stationary state, as shown in 
Fig.|2b f° r 7 = 1- For the biologically relevant case of an 
absorbing boundary, a stable stationary state does not 
exist and the cluster will always dissociate on the long 
run. In this case, one has to solve the first passage prob- 
lem of reaching the state i — for the first time. This can 
be done semi-analytically by using Laplace transforms, 
where the last backtransform has to be done numerically. 
Alternatively, one can solve the Master equation numer- 
ically by the Monte Carlo method (most efficiently with 
the Gillespie algorithm 12;] ) , as we always do in the gen- 
eral case, when both rebinding 7 and force / are finite. 
Fig-Et shows that in this case, the plateaus from Fig. 12b 
tilt downward, while po increases steadily with time t. 
Stability further decreases if force / is turned on. For 
very large force, rebinding (including the boundary type 
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FIG. 3: (a) Simulation results for the mean number of closed 
bonds N at time r for / = 0, 7 = 1 and iV = 1,2,5,10 
and 15 (lower to upper lines), (b) Four typical simulation 
trajectories for each of the cases No = 10, 100 and 1000 for 
7=1 and f/N = 0.25. Dotted lines are N(t). (c) Same 
for f/No = 0.3. (d) Comparision of stochastic (solid) and 
deterministic (dashed) results for N(t) for No = 5 and 10 for 
7 = 1 and f/N = 0.3. 

at i = 0) becomes irrelevant, because the reverse rate r 
dominates the forward rate g. Using a recursive scheme 
to construct pi from for 7 = we find 

(No \ Wo No -r(i)r 

The probability for cluster dissociation at time r is 
pi(r)r(l). Setting i = 1 in Eq. © and using Eq. J2J, 
one obtains a formula which has been given before in 
Ref. 0| . Fig. |2i shows, for the case / = 50, that now 
the cluster decays very rapidly, with only few of the in- 
termediate states being visited to an appreciable degree. 

Once the set of state probabilities Pi(r) is known, one 
can calculate any quantity of interest, in particular the 
mean number of closed bonds N as a function of time r. 
We find that N(r) usually decays exponentially. This is 
demonstrated in Fig. Eft for the case / = 0. The origin 
of the exponential decay can be understood as follows: 
first the system equilibrates into a binomial distribution 
peaked around N eq = 7^/(1 + 7) as described by the 
result for the reflecting boundary, that is Eq. JSJ. The 
lower tail of this distribution then 'leaks' into the state 
i = due to the absorbing boundary. The smaller Nq or 
7, the larger the tail contribution at i — 1 and the faster 
the systems loses realizations to the absorbing boundary. 
The resulting decay can be approximated by N(t) ks 
N eq e~ aT with a rj pi(oo) from Eq. F° r the values 
of N used in Fig. Eft, one finds a « 4.6 x 1CT 4 ,9.7 x 
1(T 3 , 0.16 and 0.5 for iV = 2, 5, 10 and 15. Numerically 
we find a = 2.5 x 10~ 4 , 8.5 x 10" 3 , 0.13 and 0.6, thus the 



leakage estimate is rather good. 

In order to assess the role of fluctuations, it is instruc- 
tive to study single simulation trajectories. Since we 
use the Gillespie algorithm for exact stochastic simula- 
tions |l2j . they are expected to resemble experimental 
trajectories. Fig. Eft shows that for large cluster size 
and small force, typical trajectories fluctuate around a 
plateau value close to N eq . However, for small cluster 
sizes, fluctuations to smaller bond numbers lead to fast 
loss of realizations to the absorbing boundary. Final de- 
cay is rather abrupt due to force-accelerated rupture for 
decreasing bond numbers. Fig. Eft shows that for suffi- 
ciently large force, also the large clusters decay quickly. 
The loss of stability for any cluster size follows from Bell's 
stability analysis of the deterministic equation Eq. 
which yields a critical force f c = N a plog(7/e) Q, where 
the product logarithm plog(a) is defined as the solution 
x of xe x — a. For typical values of No and 7, f c be- 
longs to the intermediate force regime, 1 < f < Nq. For 
7 = 1, we have f c /Na = 0.278. Fig. Eft and c are below 
and above the critical force, respectively. In contrast to 
Bell's continuum analysis, our stochastic analysis shows 
that for small clusters a small increase in loading can 
lead to the fast decay characteristic for the case with- 
out rebinding also for forces below f c . Fig. Eft compares 
N(t) as obtained from simulations to N(t) as obtained 
from numerical integration of the deterministic equation 
Eq. This shows that stochastic and deterministic 

results differ also on the level of the first moment. 

The quantity of largest practical interest is cluster life- 
time T as a function of the model parameters iVo, 7 and 
/. In general, T can be calculated from the adjoint Mas- 
ter equation [lOj. For Nq = 2 and 3, the solutions can also 
be found by directly summing with appropriate weights 
over all possible dissociation paths, each of which is a 
sequence of Poisson processes. For iVo = 2 we find 

T = i (V'/ 2 + 2e-' + 7 e- 3/ / 2 ) . (7) 

For No > 4, the direct procedure becomes intractable. 
However, in the case of vanishing rebinding (7 = 0), there 
is only one dissocation path and the exact solution is sim- 
ply T = YI}1 x l/r{i) for all values of N 0. In Fig. Hi 
we plot T as a function of / /Nq for different cluster sizes 
TVo- In the small force regime, / < 1, T plateaus at the 
value H No = J2?=i V* ~ ln7V + 1/(2JV ) + T. Here 
Hnq are the harmonic numbers and T — 0.577 is Euler's 
constant. In this regime, T depends only weakly (loga- 
rithmically) on Nq and large cluster sizes are required to 
achieve long lifetimes [T3L [l4| . In the intermediate force 
regime, 1 < / < 7V , we find T w H No - H f » ln(JV //). 
Here the effective cluster size is reduced to N /f, be- 
cause the cluster dissociates very rapidly for i < f. In 
the high force regime, f > Nq, only the term with i = Nq 
contributes: if the first bond breaks, all remaining bonds 
break within no time. The destabilizing effect of force 
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FIG. 4: Analytical results for cluster lifetime T. (a) T as 
a function of f /No for 7 = and No = 1, 10, 100, 1000 and 
10000 (lower to upper solid lines). The dotted curve is the 
approximation T = In (0.61(iVo//)). (b) T as a function of 
cluster size No for / = and 7 = 0.0, 0.1, 1.0 and 10.0 (lower 
to upper lines). 

can be counteracted by rebinding. In the case of vanish- 
ing force (/ = 0), the solution can also be found by using 
Laplace transforms [To| ■ We find 

For 7 = 0, we recover the result T = H^ a from above. 
For N = 2, we get the result T = (3 + j)/2 following 
from Eq. JJ}. In general, T scales ~ 7^0-1 w jth rebind- 
ing rate. In Fig. 0Jj we plot T as a function of TVo for 
different values of 7. For 7 < 1, the logarithmic depen- 
dence of T on Nq is valid over a wide range of cluster size. 
However, for very large clusters, lifetime starts growing 
exponentially with Nq. For 7 > 1, this strong increase 
of T with A*o is found for any value of Nq. Therefore 
increasing rebinding is much more effective than increas- 
ing cluster size in achieving cluster stability, and essential 
to ensure physiological cluster lifetimes with reasonable 
numbers of bonds. For example, in the absence of both 
force and rebinding and if the lifetime of each bond was 
1 s (ko = 1 Hz), Eq. JHJ) predicts that the astronomical 
number of ~ 10 40000 independent bonds would be needed 
to achieve a cluster lifetime of one day (T ~ 10 5 s). In 
contrast, for k on = 1 Hz (7 = 1), the same cluster life- 
time T is achieved by iVo = 20. If rebinding is ten times 
slower than unbinding (7 = 0.1), cluster lifetime T is 
down to 7 s and one needs No = 150 bonds to regain a 
cluster lifetime of one day. In this way, knowing cluster 
lifetime and two out of the three parameters Nq, 7 and 
/ allows to estimate the unknown one. 

Our model is also relevant for cell adhesion if initial 
loading is much faster than cluster lifetime (otherwise the 
assumption of constant force is not valid) and if cluster 
decay is much faster than potential reinforcement pro- 
cess (otherwise the assumption of constant cluster size is 
not valid). One example which might satisfy these condi- 
tions is L-selectin mediated leukocyte tethering in shear 



flow [l5j . Our assumptions do certainly not hold for cell- 
matrix adhesions, which have been shown to grow rather 
than to decay under the effect of force [2\. Although 
the specific processes at work at cell-matrix processes are 
not the subject of this work, our model suggests that the 
stress constant ~ 5.5 nN/ /im 2 recently measured on elas- 
tic substrates for the physiological loading of cell-matrix 
contacts through the cell's own contractile machinery [lfij 
might be close to the critical force f c — N plog(7/e), 
because then small changes in cytoskclctal loading would 
result in strongly accelerated cluster decay. Increased 
force on a subset of bonds might in turn induce signal- 
ing events (possibly through mechanical opening-up of 
protein domains) leading to subsequent recruitment of 
additional bonds. Recent single molecule experiments 
for activated as/3i-integrin binding to fibronectin gave 
k = 0.012 Hz and F b = 9 pN [l7|. Setting F c = 5.5 nN 
and using No = 10 4 , we predict 7 = 0.2, corresponding 
to a rebinding rate k on = 0.002 Hz. 
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